Influence of spin fluctuations near the Mott transition: a DMFT study 



O 



Serge Florens,^ Priyanka Mohan, ^ C. Janani,^ T. Gupta, and R. Narayanan^ 

^ Institut Neel, CNRS and Universite Joseph Fourier, B.P. 166, 
25 Avenue des Martyrs, 38042 Grenoble Cedex 9, France 
'^Department of Physics, Indian Institute of Technology, Chennai-600036, India. 
''The Institute of Mathematical Sciences, C.I.T. Campus, Chennai 600 113, India. 

(Dated: June 29, 2010) 

Dynamics of local magnetic moments near the Mott metal- insulator transition is investigated 
by a Dynamical Mean-Field Theory solution of the Hubbard model with additional fully-frustrated 
Heisenberg interactions. Appearance of deviations from pure Curie behavior in the spin susceptibility 
shows that the temperature bending of the Mott transition lines is precisely associated with the onset 
of magnetic correlations induced by spin exchange. A weakening of the effective exchange energy 
upon approaching the Mott boundary from the Heisenberg limit is also observed. 
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The Mott metal-insulator transition, wherein elec- 
tronic waves are localized by short-range electron- 
electron interactions (see [l[ for a review) , is one of the 
most complex phenomenon observed in strongly corre- 
lated electronic systems. Even though the appearance 
of a Mott gap is purely driven by the charge degrees of 
freedom, it is expected that magnetic fluctuations play a 
very crucial role in determining the true nature of this 
phase transition. In the paramagnetic Mott insulator, 
local moments are indeed well defined objects after their 
creation at high temperature (at a scale set by the local 
Coulomb interaction) and before their ultimate antifer- 
romagnetic ordering at the Neel temperature, offering a 
window in which complex behavior of the spin excita- 
tions is yet to be clearly understood. Experimentally, the 
simplest situation in this respect occurs when the low- 
temperature magnetic ordering is first order, as in the 
case of Cr-doped V2O3. Since the magnetic correlations 
are expected to be weak in this case, many predictions 
can be made from a single-site approach like the Dynami- 
cal Mean Field Theory (DMFT) [2| , where local moments 
are described as freely fluctuating from each other in the 
insulating state. As a consequence, the low-temperature 
metallic phase leads upon heating to an insulator, giving 
a Pomeranchuk-like effect where the entropy gain bene- 
fits to the state with magnetic degeneracy. Conversely, 
other classes of materials, such as the K-organics, dis- 
play a continuous magnetic transition into the Neel state, 
so that the localized spins will experience strong collec- 
tive fluctuations. The most striking experimental conse- 
quence lies in the progressive disappearance of the 
Pomeranchuk effect upon approaching the magnetic or- 
dering temperature, so that the Mott transition lines in 
the pressure-temperature phase diagram hend at low tem- 
perature. These qualitative arguments have received re- 
cent confirmation from cluster DMFT calculations of the 
phase diagram of the two-dimensional Hubbard model on 
a frustrated lattice 0-01, that take into account nearest 
neighbor magnetic exchange leading to short-range sin- 



glet formation. However, the precise connection between 
the appearance of low-energy magnetic excitations and 
deviations from the phase diagram of the single-site ap- 
proach has remained unclear. 

In order to investigate this question in a controlled 
setting, we consider the standard Hubbard model supple- 
mented with long-range fully-frustrated random Heisen- 
berg exchange interactions: 
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Note that hopping parameters t run over adjacent sites of 
a yet arbitrary lattice, while the exchange terms Jij con- 
nect all the local spin variables Si = ^, dl^^^f^d- ^, in 
a random fashion (amplitude and sign). Also, in Eq. ([T|), 
/i is the chemical potential and U is the local on-site 
repulsion, which controls the degree of electronic corre- 
lations. 

This non-trivial Hamiltonian greatly simplifies in the 
large coordination limit, where all the quantum dynamics 
can be exactly encapsulated in a single-site action: 
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Here, two different baths that depend on imaginary-time 
T have been introduced (/3 = 1/T is the inverse tem- 
perature): i) hopping processes from the local site into 
the lattice (and back) result in a retarded hybridization 
A(r — r'), which corresponds to the standard DMFT 
single site action ii) non-local magnetic coupling ex- 
actly amounts to a retarded spin interaction Q{t — r'). 
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providing friction on the local moments, hence reduc- 
ing their entropy. The latter term is easily obtained 
through a Gaussian averaging over the disordered ex- 
change coupling Jij in Eq. ([1]), followed by a dynam- 
ical decoupling of a four-spin interaction, as shown in 
S ■ For a Bethe lattice, both baths are determined 
self-consistently, through the rather simple set of equa- 
tions: 



A(t) = t^Gair) 
Q{r) = J\{t). 



(3) 
(4) 



Here, Gdir) = —{d^{T)dl,{0)) is the local d-electron 
single-particle Green's function, and xi''') = {S{t).S{0)) 
is the local d-electron spin susceptibility. The free d- 
electrons are endowed with a semi-circular local density 
of states with half-width D ~ 2t, while the degree of 
magnetic frustration is controlled by J, the width of the 
probability distribution used to perform the disorder av- 
erage. 

We now solve the still complicated quantum impurity 
problem ([2]), by laying recourse to a method that captures 
all the crucial physical ingredients. A simple possibility 
is to enlarge the spin symmetry from SU(2) to SU(A^) in 
the large N limit ll| , from which the good emergent vari- 
ables in the magnetic sector are known to be spm- 
carrying zero-charge fermions (dubbed spinous) f^, with 
a = 1 . . .N. Charge is then naturally encoded through 
a second auxiliary field, the phase </> dual to the local 
char ge, resulting in the so-called slave-rotor representa- 
tion jl3l Il4| of the physical electron: dj. = /j^e"^ = f^X. 
Following |13| . we have written the phase as a complex 
bosonic field X constrained as |X(t)P = 1 through a La- 
grange multiplier A(t), so that the action reads (particle- 
hole symmetry, i.e. /i = — [//2, is assumed here): 
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In the large N limit of [13|], the Green's functions for 
the auxiliary particles are readily obtained from Dyson 
equation in Matsubara frequency: 



G-fl{iv„) = -jj- + A - Y.x{ivn) 



(6) 

G'J^(ia;„) = «w„ - S]/(iw„), (7) 

where the associated self-energies read in imaginary time: 

Ex(r) = N/^{T)Gf{T) (8) 
E/(r) = l^{T)Gx{T)+Q{T)Gf{T). (9) 



Similarly, the parameter A in Eq. ([6|) is obtained by the 
average constraint Gx{t = 0) = 1. When re-expressed 
in terms of the emergent variables, the baths entering 
the self-consistency conditions Eqs. ((SE]) can be simpli- 
fied as: A(t) = t^Gdir) = eGf{T)Gx{T) and Q(r) = 
J'^x{t) = J^Gf{T)'^. Together with these self-consistent 
equations, the above self-energies form a non-linear set 
of integral equations that we solve numerically. We note 
that two important limits are naturally incorporated in 
this scheme: i) For ,7 = (standard Hubbard model), 
the spin bath Q(t) can be forgotten, and we recover 
the slave-rotor treatment [l^] of the usual DMFT equa- 
tions, so that a Mott metal-insulator transition occurs at 
a critical value of U/t; ii) More interestingly at [/ = oo, 
charge fluctuations are fully suppressed, as the gapped 
rotor propagator Gx{t) can be discarded in Eq. (0). 
One then recovers the quantum spin liquid equations of 
Sachdev and Ye , where the exchange J is responsible 
for changing the high temperature Curie behavior of the 
spin susceptibility xi^) — 1 into the slow temporal decay 
x{t) oc 1/|t|. As a consequence the temperature depen- 
dence of the static susceptibility x(T) = Jgdrxir) turns 
from a Curie law x(T) = 1/T at T ^ J into a slower 
logarithmic divergence xC^) ^ log(l/2^) at T ^ J. This 
interesting absence of (spin glass) magnetic ordering at 
finite temperature, due to the enhanced quantum fluctu- 
ations brought by the large N fermionic representation 
of the spin, allows us to study the onset of magnetic fluc- 
tuations in the whole range of temperatures. 

In order to investigate the interplay of Coulomb re- 
pulsion U and exchange coupling J close to the Mott 
metal-insulator transition, one needs to analyze the rela- 
tive stability of the various phases. This can be achieved 
from the knowledge of the total free energy (per site) 
which decomposes into two parts Ftot = Fimp + Fiat, 
respectively related to the local impurity problem ([2]) 
solved by the rotor equations (|8]|9]) and to the original 
lattice model ([T|) via the self-consistency equations ([Si- 
ll]). The former contribution can be obtained straightfor- 
wardly following the derivation of the slave-rotor equa- 
tion as saddle-point equations [isj : 



Finip = -A + — 2_^[lnG/(za;„) - iw„G/(iw„)] (10) 
+ ^ Y^l^nGxitiyn) - {^l/U + X)Gx{iyn)] 



-N dr[A{T)Gxir)GfiT) + -Q{T)Gf{Ty 



Taken as a functional of the pseudoparticle Green's func- 
tions, the previous expression generates the expected Ex 
and E/ self-energies Eqs. ([HE]) by differentiation with re- 
spect to Gx and G/ respectively. The lattice contribu- 
tion, associated to the DMFT equations, can similarly be 
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obtained from free energy functionals [l5|: 



Flat = 



N 



dr 



A(r)A(T) N 
i2 + ^ 



, Qir)Qir) 
dr . (11) 



Minimization of Ftot with respect to the dynamical 
fields A(t) and Q{t) reproduces indeed the DMFT self- 
consistency condition Eqs. ^SM- 

We are now equipped to investigate how magnetic cor- 
relations, controlled by the exchange term J, affect the 
Mott metal to insulator transition. To so so, we first 
construct the phase diagram in the (U,T) plane for fixed 
values of the magnetic exchange interaction J. At low 
temperatures this is done by monitoring the hysteresis 
in the local d-electron Green's function (Gdir) shows ei- 
ther slow Fermi liquid temporal decay or fast gapped be- 
havior in the metallic and insulating states respectively). 
This gives rise to two boundary lines Uci{T) and Uc2{T), 
such that a single metallic (resp. insulating) solution 
exists at a given temperature T for [/ < Uci{T) (resp. 
U > Uc2{T)). Furthermore, we determine from Eqs. (fTOl ^ 
[Tl]) the critical line Uc{T) at which the free energies of 
metallic and insulating branches coincide. The resulting 
phase diagram obtained for several values of J from J = 
to J/D = 0.5 is displayed on Fig. [T] in the case of the or- 
bital index iV ~ 3. We see that irrespective of J, there 
always exists a correlation range UciiT) < U < Uc2{T) 
bounding a coexistence region. With increasing temper- 
ature, the transition lines always merge onto a critical 
point that occurs at Tc ~ D/30 for J = 0, consistent 
with Refs. H, 13|. Interestingly, Fig. [1] shows that turn- 
ing on the Heisenberg interaction J depresses both Tc and 
f/c, leading also to a shrinking of the co-existence region. 
This can be understood by an energetic stabilization of 
the Mott phase via magnetic exchange. Turning to the 
zero temperature limit, one recovers by the collapse of the 
Uc{T) and Uc2{T) lines at J = the continuous metal- 
insulator transition described by the Brinkman-Rice sce- 
nario [H, [l3l ■ In that case, the metal is the true ground 
state for all U < Uc2{T = 0). However, as seen in Fig. [U 
the switching of J markedly separates the actual criti- 
cal line Uc{T) from the two stability lines, so that the 
Mott transition becomes first order at zero temperature, 
in agreement with recent findings Q. 

Next, we turn our attention to the orientation of the 
transition lines as the interaction between the local- 
moments is increased. At J = 0, the Uc{T) line is such 
that heating the correlated metal always leads to the 
Mott insulator, a Pomeranchuk-type effect that relates 
to the high entropy of the insulating phase. This feature 
persists even above the transition temperature: by deter- 
mining at fixed T > Tc the crossover interaction Uc{T) 
at which the d-electron Green's function (taken e.g. at 
T = 13/ 2) displays an infiection point, one distinguishes 
between a bad metal (i.e. an increasing resistivity with 
increasing temperature with non quadratic behavior) and 
a bad insulator (i.e. a non-gapped decreasing resistivity 
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FIG. 1: Metal-insulator transition phase diagram at A'^ = 
3 for several values of J = 0,0.1,0.2,0.4,0.5 (right to left). 
Labels of the various critical lines, as defined in the text, are 
explicited for J = 0. 



with increasing temperature), see pi5| . The continuity of 
the critical Uc{T) line (determined from the free energy 
at T < Tc) into the crossover line Uc{T) (determined from 
the Green's function at T > Tc) serves as a vindication 
of our methodology. We also note that such a crossover 
at T > Tc is seen experimentally in the magnetic and 
transport measurement of k-BEDT salts [J. Now, as 
J is increased one expects the entropy of the paramag- 
netic Mott state to be reduced by magnetic exchange, 
thus weakening the Pomeranchuk lineshape. This is in- 
deed seen in Fig. [1] by a progressive leftward bending of 
the critical lines at increasing J, and consistent with re- 
sults of Ref. 0, B A particularity of our model is the 
re-entrance of the Pomeranchuk effect at the lowest ob- 
servable temperature, which we attribute to the reduced 
(with respect to log(2) of free moments) yet non-zero 
entropy of the insulating spin liquid state |lol . [l2| . We 
have checked also that increasing orbital degeneracy with 
larger N values tends to weaken all magnetic effects, re- 
covering the standard J = phase diagram. 

A key point of our study concerns the connection be- 
tween the exchange induced modifications of the phase 
diagram (which can be observed e.g. in transport exper- 
iments) and the actual spin dynamics. The spin suscep- 
tibility xi''') in the absence of exchange J = reduces 
trivially to a pure Curie law in the whole Mott state 
U > Uc2, as checked in Fig. [51 At finite J, Gurie behavior 
is seen only at high temperatures such that T ^ J, while 
departure from the independent local moment picture is 
observed by a slower logarithmic divergence at zero tem- 
perature. A crossover scale Tsi characterizing the onset 
of the spin liquid state can be accurately determined from 
a fit of x{T) by the appealing form (l/Tsi) ln(l + T^i/T) 
which follows the expected high and low temperature 
regimes (the agreement of this Ansatz with our data is ex- 
cellent). This procedure allows us to extract, for a given 
value of J, the {/-dependence of the spin-spin correla- 
tion scale Tsi , see Fig. [31 obtained for = 3 and several 
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FIG. 2: Spin susceptibility xiT) versus temperature for J = Q 
(circles) and J = 0.2 (squares) at A'^ = 3 and large U . A 
l/T Curie behavior applies to the whole J = curve (free 
moments), while a spin-liquid Ansatz (1/T'3;)ln(l + Tsi/T) 
describes the case of finite J — 0.2 exchange (full line). 



J values. Wc note first that Tgi saturates at very large 





U/D 



FIG. 3: Superposed plot of the [/-dependent spin- liquid on- 
set temperature Tsi for several J = 0.1,0.15,0.2 (circles, bot- 
tom to top) with the corresponding metal-insulator phase di- 
agrams (right to left). 

U values to (half) the bare magnetic exchange constant 
J, an expected result since only spin dynamics plays a 
role upon freezing of the charge excitations, leaving J as 
the only relevant scale. More interesting is the net de- 
crease of Tsi when U is diminished towards the critical 
Mott boundary. This feature signals a weakening of the 
spin correlations in the vicinity of the Mott transition, 
brought by the finite probability of double and zero occu- 
pation of the d state (valence fluctuations in the particle- 
hole symmetric regime) . Quite intriguing is also that the 
spin susceptibility continues to display the characteristic 
signature of the spin liquid even in a regime where the 
ground state is metallic, close to the coexistence region. 
This result is reminiscent of a previous study by Par- 
collet and Georges |12|, where the doping- induced Mott 
transition (at U — oo) within the fully random Heisen- 
berg model was considered. These authors showed that 
the magnetic crossover scale Tgi persists at small doping, 
leading to a slush state of incoherent metal with spin liq- 
uid correlations, similar to what we observe here at half- 
filling in the intermediate Coulomb interaction regime. 



Our last yet most important prediction is the precise 
connection of the spin liquid scale T^i to the tempera- 
ture where the bending of the metal-insulator crossover 
takes place, as manifest from Fig.[3l This comparison in- 
deed shows that magnetic exchange is responsible for the 
quenching of the Pomeranchuk effect and the bending of 
the Mott transition lines. 

In conclusion, we have investigated the role played by 
fluctuating local moments coupled via frustrated long- 
range interactions near the Mott transition. The metal- 
insulator transition was found to be modified by the 
Heisenberg exchange in favor of the insulating state and 
even to turn discontinuous at zero temperature. Study- 
ing the magnetic susceptibility, we were able to connect 
the bending of the critical lines to the appearance of col- 
lective spin correlations, with reduced effective exchange. 
These effects should be generic to other microscopic mod- 
els and experimental systems near the Mott transition. 

We would like to acknowledge helpful discussions with 
S. R. Hassan, R. Moessner, and T. Vojta. 
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